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Abstract 

Motivated by recent studies on the dynamics of colloidal solutions in narrow channels, we consider 
the steady state properties of an assembly of non interacting particles subject to the action of a 
traveling potential moving at a constant speed while the solvent is modeled by a heat bath at rest 
in the laboratory frame. Since the description, we propose here, takes into account the inertia 
of the colloidal particles it is necessary to consider the evolution of both positions and momenta 
and study the governing equation for the one-particle phase-space distribution. We first derive 
the asymptotic form of its solutions as an expansion in Hermite polynomials and their generic 
properties, such as the force and energy balance and then we particularize our study to the case 
of an inverted parabolic potential barrier. We obtain numerically the steady state density and 
temperature profile and show that the expansion is rapidly convergent for large values of the 
friction constant and small drifting velocities. The present results on the one hand confirm the 
previous studies based on the dynamic density functional theory (DDFT) when the friction constant 
is large, on the other hand display effects such as the presence of a wake behind the barrier and a 
strong inhomogeneity in the temperature field which are beyond the DDFT description. 

PACS numbers: 82,70.Dd,61.20.-p,05.70.Ln 
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I. SHORT INTRODUCTION 



In recent years we have witnessed the emergence of a new branch of applied physics 
named microfluidics, which is the science of designing, manufacturing devices and processes 
that deal with volumes of fluid on the order of nanolitersi^^. Microfluidic systems have 
diverse and widespread potential applications^^. Some examples of systems and processes 
that might employ this technology include ink-jet printers, blood-cell-separation equipment, 
biochemical assays, chemical synthesis, genetic analysis, drug screening, electrochromatog- 
raphy, surface micro-machining, laser ablation, and mechanical micro-milling. Not surpris- 
ingly, the medical industry has shown keen interest in microfluidics technology. 

Such advances in manipulating fluids*^*!! have recently motivated Penna and Tarazona^ 
to consider a model representing a simple device to push a dilute solution of colloidal particles 
along a narrow channel. In particular they studied the effect of a moving barrier on a 
system of non interacting colloidal particles described by overdamped Langevin dynamics. 
Under the action of the potential barrier shifting at a constant speed, the fluid achieves a 
steady state, with density distribution and local current following the moving barrier. These 
authors showed that such a steady state can be conveniently studied within the DDFT— Qddk 
formalism, since the structure of the relevant equations becomes similar to that of the Euler- 
Lagrange equations describing a fluid at thermodynamic equilibrium^. 

On the other hand, the present authors in a recent paper—, hereafter referred as Ref. 
I , have considered how the inertia of the particles may modify the DDFT picture. They 
assumed that the colloidal particles have inertia, i.e. are governed by a second order stochas- 
tic equation. The governing equation for the associated phase-space distribution turns out 
to be the Kramers equation^ and represents the evolution of both positions and momenta 
of the particles. Since such a representation is still too complex and often redundant, the 
authors considered a contraction of such a description by rewriting the Kramers equation in 
terms of the infinite hierarchy of equations for the velocity moments of the phase-space dis- 
tribution. In ref. I the hierarchy was truncated systematically by means of a multiple time 
scale technique, which lead to a self-consistent equation involving only the one-body density. 
This equation is similar to to the DDFT equation, but contains additional terms taking into 
account the presence of momentum and energy currents. While in ref. I we considered only 
transient effects, namely the decay of initial perturbations towards the equilibrium, time 
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independent, state, in the present work we illustrate how the inertial dynamics affects the 
behavior of systems in situations in which a steady state is induced by the presence of an 
external time-dependent potential. The results of the present paper show pronounced dif- 
ferences with respect to the DDFT study of Penna and Tarazona and in particular display 
an exponential decay in the structure of the density profile behind the barrier which was not 
predicted by the DDFT. Moreover, we find that also the local temperature is non uniform 
throughout the system due to the heating produced by the barrier. 

We believe that these findings are generic to non equilibrium systems where the equilibra- 
tion mechanism provided by the heat bath is not very rapid. We have shown that when the 
friction is not sufficiently high the density alone is not sufficient to characterize the steady 
state of the system, and additional fields are necessary to provide a complete description. 

More generally speaking, we believe that the use of the DDFT, is justified when the 
currents are of diffusive character, while in the cases where convective terms are present it is 
necessary to include extra terms which describe the transport of momentum and energy^ 1 ^ 

The present paper is organized as follows. In section [Til after presenting the model we 
give the structure of the general solution of the Kramers equation in the region where the 
potential is vanishing small. In IHH we specialize the treatment to the case of steady state 
conditions and derive explicitly the behavior of the phase space distribution in the region 
where the traveling potential vanishes. We also derive the relation between the total force 
exerted by the barrier on the particles and the friction due to the bath. Finally in IIVI we 
give explicit numerical solutions of the Kramers equation in the case of an inverted parabolic 
barrier. We conclude the paper with a short discussion in [V] 

II. KRAMERS EQUATION FOR SHIFTING POTENTIAL BARRIERS AND ITS 
FREE MODES 

The problem of the steady states in a fluid, under the action of a shifting exter- 
nal potential, has been considered within the DDF under several conditions and model 
interactions^ 1 ^ 1 ^. In all these treatments the inertia of the particles did not play any role. 
Here we wish to consider how the inertial effects modify that picture, and to such a pur- 
pose we consider here the simplest case, which could describe a dilute solution of colloidal 
particles dragged along a narrow channel under the action of a moving potential barrier, 
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modeled by a time dependent external potential, V ext (x, t) = V ext (x — ct), which acts on the 
colloidal particles but has negligible effects on the solvent. To such purpose we consider here 
an assembly of non interacting identical particles of mass m moving in one dimension, and 
described by the following stochastic dynamics^^ 

d?x dx 

with a bath providing the particles a friction constant 7, and a thermalizing noise with 

(C(m(s))^2 im k B T S(t-s), (2) 

at temperature T Q . The external force associated with the traveling potential is f ex t{x,t) = 
— 4-V ext {x — ct), and the properties of the system can be studied by considering the equa- 
tion governing p(x,v,t), the density distribution in phase space of a single particle. The 
associated Kramers equation^ 1 ^ reads, 



d 



dt 



d , fext{x - ct) d 



—p(x,v,t)+ v— — I — p(x,v,t) = 'y ^-H jr^ p(x,v,t) (3) 



d T 8 



,2 n 



dv m dv 2 



We assume that the shifting external potential is localized within a finite region and van- 
ishes outside. Therefore, far away from such a region we should have a time independent 
equilibrium distribution p (x, v) = p exp(—v 2 /(2v^))/(\/2TTVT), where p a is the density of 



particles and vt = a/ ksTo/m the Gaussian width for their velocity distribution. For a static 
external potential, i.e. the c = limit of flTJ), the distribution p(x,v,t) would evolve in time 
towards the thermal equilibrium value p eq (x, v) = p (x,v)exp(—V ex (x)/kBT ), which would 
be reached (sooner or later) from any initial distribution p(x, v, 0). For the continuous 
shift of the external potential implies a permanent perturbation of the thermal equilibrium, 
but still there would be a transient evolution from any p(x, v, 0) to a unique stationary state 
p(x — ct, v) in which the time dependence is reduced to a shift of the x coordinate, to follow 
the external potential V ext (x — ct). This steady state is the object of the present study. All 
the results presented here may be translated to a purely static distribution in the presence of 
a time independent external potential V(x) = —f x + V ext (x), with a constant slope plus the 
same potential barrier which we are considering in eq. Q. The time derivative in the first 
term of eq. ([3]) vanishes, but there is an extra term proportional to f Q , to take into account 
the constant background force added to the localized barrier force f ex t{x). Away from the 
barrier the particles move at constant mean velocity v Q = f / (rwy), and a change of reference 



framework from v to v' = v — t> , leads to exactly the same equation ([3]) for p(x, i/, t), when 
the barrier appears as moving at rate c = —v D . These two equivalent versions of the same 
problem have been studied within the DDF formalism^ 1 ^ 1 ^ 1 ^, valid for large 7. The same 
exact mapping between the moving barrier in a flat background and the static barrier in a 
sloped potential would be valid for partially damped systems explored here. 
It is convenient to introduce the following dimensionless variables: 



T = tV T po, V EE 



Vt 



X EE X p 



c = 



Vt 



(4) 



r = 



7 



F ext (X,r) = fext{X 2 Ct \ P(X,V,r) = ^p(x,v,t), (5) 

Po 



v T Po mv T p 

Accordingly, Kramers' evolution equation for the phase space distribution function can 

be rewritten with the help of relations (jUl5]) as: 



l aP(yr) = LfpP(x? Vj t) _ Ap (X , V, r) - ±F e (X, r)-^P(X, V, r) (6) 

having introduced the "Fokker-Planck" operator Lpp whose eigenfunctions H^(V) have the 
property: 



L FP H^V) = 



_d_ 
W 



_d_ 

W 



+ V 



H,{V) = -pH,{V), 



for p = 0, 1, ...,and have the explicitly representation: 

1 ft? 1 

HJV) = -^(-lY^—exp(--V 2 ). 



(7) 



(8) 



It is convenient to define raising and lowering operators in the eigenfunctions series, 
a±H l j / (y) = i7 /i ±i(V r ), so that the contributions of the damping and the external forces in 
the last two terms of eq. ([6]) may be represented through 



VH^V) = H^ +1 + pH^(V) = (a+ + pa^H^V), 



(9) 



and 



dV 



-H, +1 (V) ee -a+H^V) 



(10) 



The exact solutions of eq. ([6]), in the regions where the external force vanishes may be 
written in terms of the infinite series of modes, p = 0,1, with the generic forrn^ 

a_L d 



P M (X, V, r) = exp(-/TY) exp 



r dx 



1 + rScY H ^ v)(t)M{x ' T) 
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The function 0^(X, r), which fully defines the mode P^'(X, V, r) represents any solution 
of the diffusion equation 

^-(p M (X, t) = -^—(h^(X, r). (12) 
From (jlip in the case /x = we obtain explicitly 

P<°>(X, V, r) = flbCVO^C^, - ^f^l) + + .... (13 ) 

which describes a density inhomogeneity, represented by the term 4>°(X, r), and the associ- 
ated momentum current, the term of order 1/r, kinetic energy current, the term of order 
1/r 2 , and so on. These terms are slaved by the density and their shapes are given by the 
successive derivatives of <fi°(X, r) with respect to X. Similarly, from ffTTl) the solution with 
/i = 1 has the explicit representation 



P (1) (X,V,t) = exp(-rY) 



H 3 (V)d^(X,r) 

2!r 2 dx 2 
+ f r o(y) ^x f 9x^ + - 



where the first line in the r.h.s. has the interpretation of a master current inhomogene- 
ity r), which slaves higher order moments with decreasing amplitudes (1/r,...), 
while the second line in the r.h.s. has the same structure as PW(X, V,t) with amplitude 
0(o) _ Y~ 1 dx4>^ 1 \ an d both terms have the fast decay of the exponential pre-factor. The 
physical interpretation of such a combination is that an initially pure current fluctuation, 
described by Hi(V)4>i(X, 0) would die very fast, as exp(— Tr), but leaving behind a density 
fluctuation proportional to T~ 1 dx4>^ 1 ' (X,0), which would evolve diffusively. The particular 
combination in (j!4p is such that it completely cancels that remnant density fluctuations, 
i.e. it orthogonalizes P^(X, V, r) to P(°'(X, V, r), and leaves a purely fast decaying form. 
The generic free mode of order /x, is a master term ^»w(X, r)H^(V)exp(— IV), representing 
a density (/x = 0), current (/x = 1), temperature (/x = 2), heat (/x = 3), etc., perturbation 
of the equilibrium distribution p Q (x,v). The master distribution (p^)(X, r), slaves the per- 
turbation components associated to any other H^iV), with increasing powers of the inverse 
damping 1/r, so that whole distribution P^>(X, V, r) decays towards equilibrium with an 
exponential decay time (/xr) _1 . For time independent external potentials, the high order 
modes are only visible as very short transient states of P(X,V,t) towards p (x,v), and 
the in the large damping limit, T ^> 1, the modes are essentially reduced to their master 



component^. We analyze in this work the role of these modes under the continuous shift of 
the external potential, for finite values of the damping constant T. 

III. STEADY STATE SOLUTION 

A. The steady state form of the free modes 

Now, we impose the steady state condition P(X, V, r) = P(X — Ct, V), shifting with time 
to follow the boundary conditions in the moving potential barrier, in terms of the variable 
X = X — Cr.i We analyze first the form of the free modes of the expansion to represent the 
solution of eq. ([6]) in the regions where the external force vanishes. Since the steady solution 
has the property 

^ [PM(X,V,T)] = -cA [pM(x, V»] , (15) 

it follows that 

|- [exp(- M rr)0^(X,r)] = -C exp(-AiIY) A [^(X,r)] , (16) 

so that we can transform the diffusion equation ([T21 for the master distribution into an 
ordinary differential equation 

(>1 4>W{X,t) +rcA M(x,r) -//r 2 0W(X,r) = 0, (17) 



dX 2 ^ v ' ' ox 

whose solutions are proportional to exp((3± X) with 

_-C± VC 2 + 4/i 



(18) 



r 2 

Finally the product exp (-/ilY) <f>M(X,r) featuring in eq. (TlT!) has the form consistent 
with eq. ( fTBT) : 

exp(-^rr) <p M {X, r) = J2 4^ ex P l^i* ~ C T )\ ' ( 19 ) 

The amplitudes A±^ determine the contribution of each mode in any region where the 
external potential vanishes. Since we assumed that the potential barrier is restricted to a 
finite region around X = X — Ct w 0, we shall refer as the front region to the positive 
values of X, whereas negative values X represent the wake region. 
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For the first mode, p = 0, the exponent (3 + vanishes, so that A + = 1, to represent the 
only possible constant contribution to P(X,V,r), the equilibrium distribution p(x,v,t) = 
P H (v/vt)/vt, away from the perturbation. The second exponent for /x = 0, is (3^ = — TC, 
so that it can only contribute to P(X — Cr, V) in the front region of the advancing potential 
barrier, with an amplitude A_ to be fixed by the boundary condition at the advancing 
front of the external barrier. On the region left behind the barrier, the amplitude A_ 
has to vanish, since otherwise P(°'(X — Cr,V) would diverge as exp(— CX) for X C 0. 
Therefore, substituting the solution ffT9l with /i = into eq. f|T3|) we obtain the structure 



P(°\X,V) = H (V) + A { °>e 



(o) -rex 



C 2 H 9 (V) 
H (V) + CH 1 (V) + — ^r J - + - 



(20) 



for X = X — Cr on the front side of the advancing barrier, while behind the barrier we have 
the pure equilibrium structure p(°\X, V) = Hq(V), with no remnant wake structure. 

The contribution proportional to H (V) in P^{X,V) has precisely the shape obtained 
from the analysis of eq.([T]) in the strong damping limit 12 , when the particles are always at 
their limit velocity and the inertial term can be neglected. In this limit the Smoluchowski^ 
description of the system is sufficient, and the solution can be written as P(X,V,r) = 
p(X,t)H (V)/p , where p(X, r) satisfies the following diffusion equation with drift 

dp ^ T) = - i A ( P (X, r)F e {X - Cr)) , (21) 

dr T dX 2 r dX y ' 

and the stationary solution p(X — Cr) for shifting potential barriers, has the exponential 
front and the complete lack of wake identical to the Hq(V) contribution to ([201 . The 
only qualitative difference between the fully damped system described by the Smoluchowski 
equation, and the \i = mode solution, of ([2D1 is that the front density perturbation slaves 
a current CHi(V), a kinetic energy increase C 2 Fi2(V)/2, and similar higher order terms 
which may be resumed to give exactly the form 

P (0) (X, V) = H (V) + A^ ] e- Tcjc H (V - C), (22) 

i.e. the whole perturbation oi P^(X, V) over the equilibrium value H (V) has a Maxwellian 
distribution of velocities but shifted to the reference frame of the advancing potential barrier. 

All higher order terms are characterized by (3^ > and [3^ < 0, so that the exponent 
j3^\ has to be taken at the front side and (3^ behind the barrier, so that there is one free 
amplitude A±^ for each mode at each side of the barrier. The distribution functions for 
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these modes may also be written in terms of the shifted eigenf unctions of the operator Lpp, 
H V {V + /3± /r). Thus, the /x — 1 mode has the form 



whereas for the /x = 2 mode we find: 



H 1 [V + 



r 



+ 



Ho V + 



r 



(23) 



P®(X,V) = A%>e$ > * 



and the generic structure of the /x mode is 



V) = ( 1 + " hJv + ^ 



(24) 



(25) 



Notice that all contributions P^(X, V) for /x > decay exponentially away from the barrier. 

The inclusion of higher order terms creates a wake density fluctuation structure, with 
exponential decays exp(/3+ X), which have /3+ ~ y/Ji for C <C 1 and /3+ ~ /x/C <C /x for 
C ^> /x. The front density structure contains several exponential decays exp((3_ X), with 
/3_/r ~ — ^//x for C < 1 and (3+^ /Y ~ — C for C ^> /x. Both at the front and the wake 
regions, the density fluctuations go together with fluctuations in the velocity distribution, 
which may be described as shifted equilibrium distributions, H (V + /3± ), shifted current 
distributions Hi(V + /3± ), etc... The front region is broad if the damping is weak and 
the barrier velocity small, because the restoring force is proportional to the velocity of the 
colloidal particles with respect to the quiescent solvent. The velocity distribution changes 
in front of the barrier and develops secondary peaks at V = C, V = —j3_ /Y : V = — (5_ /Y, 
etc. 



B. Generic properties of the steady state produced by a shifting barrier 

We consider now the generic solution of eq. ([6]), including the regions inside the moving 
barrier, where we have to include the force term. The steady state condition P(X, V, r) = 
P(X - Ct, V) = P{X, V) transforms eq. © into 

{v - o s -mn = rL FrP( x, v) - (26) 
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The general solution of this equation may be represented as 

oo 

P(X,V) = Y t $ v (X)H v (V), (27) 

with generic functions Q U (X), to be determined from the projections of eq. (1261) on each of 
the FP eigenfunctions H U (V). The projections for v — and v = 1 give 

M-c^E! (28) 

and 

dx dx dx 

The general form for any v > 1 is 

( „ + d«W*) _ + ^« = - „«.(*). (30, 

aX aX oX 

In absence of the force term F e (X), the general solution of this (infinite) set of coupled 
ordinary linear differential equations may be written in terms of the free modes ( 125|) . with 
arbitrary amplitudes at the back side, and at the front side of the moving barrier. 

The structure of eq. (1281) is independent of F e (X), and it represents the continuity 
equation, relating the mass density p(X) = p o $ (X) to the current density j(X) = p a <&x(X), 
to keep the mass balance under a steady flow, 

dj(X,r) _^ dp(X-Cr) 

~^x~~ c dx • 

The integration of (1251) from the boundary conditions <&q(X) = 1 and $i(X) = 0, far 
away from the moving barrier, gives 

$ 1 (X)=C($o(X)-l) ) (32) 

i.e. any positive excess <&q(X) — 1 > in the distribution of particles near the moving barrier 
is associated to a current j(X) = C{p(X) — p ) following the barrier shift. The regions with 
®a(X) < 1 imply a depletion of the density, and a counter-current with opposite sign to the 
barrier displacement. In the strong damping limit 12 such a depletion and counter-current 
were limited to the interior of the potential barrier, since there was no wake left behind 
it. The inertial effects here included open the possibility of such wake, so that we may 
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find regions outside of the moving barrier where the mean velocity (V) = <3>i(X)/$o(X) = 
C[l — l/$o(X)] has the sign opposite to C. 

Eq. ( 1291) . from the projection of eq. (1261) on H\(y), represents the local balance of 
momentum. If we integrate it across the whole inhomogeneity, from far from the rear to 
far from the front of the moving potential barrier, the integrals of all the derivatives vanish, 
and we get that the total force F T , produced by the barrier on the particles balances the 
friction force created by the bath on the total current 

/OO P "X 

dXF e (X)%(X) = T / dX^iX), (33) 
-oo J — oo 

i.e. it gives the global force balance in the system. Notice that only the region of poten- 
tial barrier contributes to the first integral in the left hand side, while the entire volume 
contributes to the right hand side. 

The integration of ([29]) from <3> P0 = 1, $iPO = 0, and $ 2 (X) = 0, at any point far 
from the barrier gives the local excess of kinetic energy at any point, 

-x 



C$x(X) - $ (X) - 1 + ! dX (f £ (X)$ (X) - r$ x (x) 

J —oo 



(34) 



Therefore, once we have the particle distribution $o(X), we may get the mean velocity of 
the particles (V) = <3>i(X)/$ (^0 from fl3"2l . and their local temperature relative to that of 
the bath, T(X)/T = 1 + $ 2 (X)/$ (X) from $M§. 

Similarly, the equation for v = 2 in the series (|30|) corresponds to the energy balance. Its 
integration from a point far behind barrier to an arbitrary point X gives direct access to the 
heat current $3(X), while its integral across the whole inhomogeneity gives the total power 
transferred from the barrier to the particles 

/oo poo 
dXF e {X)^ 1 {X) = 2T dX$ 2 (X), (35) 
-oo J —oo 

where the last integral has to be interpreted as the total heat dissipated by the particles 
due to the local temperature difference over the bath, $ 2 = $o(X)(T(X)/T — 1). Notice 
that the steady state conditions, and the fact that the potential energy vanishes both at the 
front and at the rear of the moving barrier, gives a direct relation, W = OF?, between the 
power and the force. Through eqs. ( J32l) . ( J33l) and ( !35l) we get also a relationship between 
the total excess of particles and the excess kinetic energy. Written in terms of the original 
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variables, 

dxp(x)(T(x) - T ) = — / dx(p(x) - p ), (36) 

this should be a generic property of the steady state distributions, independent of the damp- 
ing r. 

C. Expansion in terms of the steady free modes 

The above expressions for $i(Jf), $ 2 pf),.., given in term of <&q(X) can only be used 
after the whole set of ordinary differential eqs. (|28|) - (|30|) are solved. That requires either 
a resumation of all the terms, as done for the free modes in eq. (1221) . or some truncation 
scheme to perform a numerical integration for the regions with F e (X) 7^ 0. Unless the force 
is very weak everywhere, a direct truncation scheme of the expansion in eq. ( 1271) . e.g. taking 
$ 3 (X) = 0, and solving the first three equations to get <& (X), &i(X), and $ 2 (X), leads 
to unphysical results, strongly dependent on the order of the truncation. On the contrary, 
we have found very good convergence, at least for any T > 1, using a finite parametrization 
of P(X, V) based on the natural modes for the free particles. We fix the number p ma x of 
such modes to be used in the front and in the wake regions, so that the solution, P(X, V), 
is described by p ma x + 1 constants at the first region, and p ma x constants at the 
second region, besides the trivial contribution Ay = 1. Within the barrier region we use 
2/i m aa; + 1 independent functions, ip±(X) to parametrize P(X, V) as 

H*,v) = g £ ^ ^) h, (v+^)- (37) 

Therefore, each term in the expansion (j27|) is expressed as a linear combination of the 
functions ijj±\x) to be determined by means of eqs. (1281) - (I30~]) .for all values v < 2p max + 1. 

The simplest parametrization within this scheme corresponds to include only the /i = 
mode, with (3^ = and f3_ = — CT, so that 

P(X, V) = ^(X)H (V) + ^\X)H (V - C). (38) 

Hence, all the terms in expansion f l27j) are given in terms of these two functions, 

$ (X) = ^\X)+^(X), = C^\X), $ 2 (X) = ^\X), etc., 

(39) 
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The projections of eq. (I26I) on the first two FP eigenf unctions are enough to determine 
*IJ { +\X) and ^{X). From eq. (1281) we get that ip+\x) has to be constant all over the 
system, both inside and outside the potential barrier, therefore it is fixed by the asymptotic 
value if)+\x) = 1, and we may use the particle distribution $> (X) = 1 + tp_\x) as the 
only free functional variable. Regarding now the projection on Hi(V), we get that the 
contributions from the derivatives of <&i(X) and $2p0 on the left hand side of eq. (|29|) 
cancel each other, so that 

= F e (X)%(X) - CT(%(X) - 1). (40) 

oX 

which is exactly the DDF equation obtained and solved by Penna and Tarazona^ from the 
integration of (I2T]) . Notice that this simplest parametric description of P(X, V) is therefore 
consistent with respect to the mass and momentum balances, but it has not the flexibility 
to recover the equivalent balances of energy [y = 2), heat current (u = 3), etc., required 
by eqs. ( 1301) . A direct substitution of eqs. ( 1391) into eqs. ( 1301) shows that the local balance 
for v > 2 fails by a term C v ~ x F e (X) /{y — 1)!, at each v > 2. Such a failure is less important 
for low shifting rate, C< 1, and for modes v ^> C. Also, the global balance represented by 
eqs. fl33|) - fl35|) would be kept at any order u, since the total integral of F e (X) has to vanish. 
The inertial effects appear to recover the local balances missed by the DDF approximation, 
and we may include them in a systematic way including in P(X, V) the contributions of 
the higher order free modes. That enlarges the set of free functions i()±(X), and allow the 
solution of eq. (|26|) up to higher order eigenf unctions of the FP operator. 



IV. NUMERICAL RESULTS FOR A PARABOLIC POTENTIAL BARRIER 

As an application we study a parabolic potential barrier of the form, U(X) = k(1 — X 2 )/2 
creating a linear force F e (X) = kX, restricted to the interval — 1 < X < 1. As we consider 
only the steady case we have to solve eqs. ( !30l within the barrier, and to find the solutions 
matching with the physical solutions eq.( l25l) at the front (X > 1) and at the wake (X < 
— 1). The matching of P(X,V) inside and outside the barrier is achieved by requiring 
that P(X, V) in eq. fl26|) has to be continuous at X = ±1, but with a discontinuous first 
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derivative with respect to X, to match the discontinuity in F e (X), i.e.: 



(V-C) 



( dP(X } V) 



dX 



dP(X, V) 



dP(X, V) 

K 



X=l-e/ 



dV 



(41) 



x=i 



and a similar condition at X = — 1. Most of the results presented here have been obtained 
with Umax = 4, i.e. with nine independent functions ip±(X), besides the trivial ip^ — 1j 
to ensure the correct projection of eq. (125]) up to order Hg(V). Nearly identical results are 
obtained with \i m ax = 3, and even with /i mai = 2 for T > 2. However, qualitative differences 
appear with respect to the DDF result (p, max = 0), unless we have both large damping V 
and a low shifting rate C for the barrier. In any case we have to deal with a set of linear 
differential equations, with a shooting boundary problem, to get the physical match with 
the free modes, so that A { f = at X = 1, and A M = at X = -1 

In Fig. |V] we present results for a high potential barrier, k = 10, moving with respect 
to the bath at a relatively low velocity, C = 0.2. For large damping the system is in the 
strong drift limit 12 . The density distribution is strongly depleted within the barrier, while 
the density in the front region grows to a large value, more than sixty times the asymptotic 
density in this case, so that there are enough particles going over the barrier to keep the 
stationary state. We observe that for V > 1, the inertial effects have little influence in the 
structure of the front. When $>q(X) is rescaled in terms of Y(X — 1), as in Fig. IVTa). 
the curves collapse into a single large T limit. This is consistent with the fact, that the 
velocity distribution at the front region is dominated by the shifted Maxwellian form fl38l) . 
The effect of reducing T below the value 1 renders smaller the amplitude of the exponential 
contribution in the formula <&q(X) = 1+Aexp(—CTX). Nevertheless, for the lowest value of 
T presented in that figure the expansion in modes is still far from convergence for /i max = 4. 

In Fig|V](b) we present the structure of the wake by rescaling the distance from the left 
edge of the barrier by the the factor YC. The profile saturates for low Y, while is continuously 
reduced as Y increases. This is consistent with the no-wake prediction in the large damping 
limit, when the inertial effects are fully suppressed. Nevertheless, the decrease of the wake 
structure with increasing Y is very slow, so that the presence of such region, with p(x) < p 
and hence mean velocity (v) = vt(1 — p / p{x)) < 0, behind the shifting potential barrier, is 
an important qualitative effect induced by the inertial dynamics of the particles, and which 
was neglected within the DDF analysis^. 

In Fig. |V]we present the results for the same barrier as in Fig. |V]but with a much larger 
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velocity C = 2. In the large damping limit such a situation corresponds to a "high counter- 
current" regime^, in which the barrier moves too fast to produce a strong perturbation 
in the density distribution. When the bath damping parameter is reduced, the inertial 
dynamics creates a strong amplification of the front structure, which is now much more 
symmetric with respect to the advancing barrier front at X = 1. Roughly speaking, a half 
of the particles at the front are actually within the potential barrier, for ^ X < 1. This 
is to be compared with the result at low C where the advancing front was mainly located 
at X > 1. The density depletion is limited to the rear edge of the barrier X pa — 1, and the 
density is never lower than 0.5 the asymptotic value. The structure of the wake region does 
not show the strong scaling effect with T observed in Fig. |V] (b), for the low C case. The 
scaled structure of the front in terms of TC(X — 1) is presented in Fig. IVTa). and shows that 
the decay of the density is still well represented by the exponential form A_ exp(— TCX) of 
the zeroth order mode, but with a T-dependent amplitude A_. The maximum amplitude 
of the wake, just behind the moving barrier seems to be similar for all the cases with large 
T, while the decay increases with T. The results in terms of the scaled distance TC(X + 1) 
may be compared with those in Fig. |V](b) for the slowly moving barrier, and we observe 
that the wake extends now further away from the barrier edge. 

We turn, now, to the study of the local rescaled temperature^ obtained from eq (1511) . 
The results for the C = 0.2 case in Fig. |V] indicate that the shifting barrier produces a very 
strong heating of the system within the barrier, with maximum T(X) pa 25T D at the rear 
side of the barrier, in the region of lower density. At the scale of the maximum T(X) the 
temperature is apparently constant at the front side, but the inset shows that there is a 
sharp rise of temperature at X < 1, and also we observe a kind of precursor plateau, over 
distances of the order X pa 20/T from the barrier edge, and with a T independent value 
T(X)/T pa 1.5. The width of that plateau may be understood from the huge enhancement 
of the density on the front side of the barrier, so that until the exponential decay makes 
&o(l)exp(— TC(X — 1)) pa 1, the large majority of the particles contributing to <&o(X) 
belongs to the exponential component of that front, and the value of T(X) at the plateau 
would represent the temperature of the advancing front. The structure of T(X) the wake 
is much narrower than at the front, and it indicates a moderate heating within the density 
depletion shown in Fig. |Vj 

In Fig. |V] we present the temperature distribution for the high velocity case, C — 2, 
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described in Figs. |V] and |V] The maximum temperature is T(X) m 5.5T D , and it is still 
located at the rear half of the barrier, associated to the minimum density. The "precursor 
film" at the front is much shorter and higher, so that only for the lowest value of T may may 
be interpreted as a incipient "plateau", this is consistent whit the interpretation given above 
when we consider the density distributions in Fig. |Vj The most peculiar feature of T(X) 
at this high value of C is the appearance of local minimum between the main maximum 
and the front edge of the barrier. The relative importance of this feature increases with 
decreasing T, i.e. as the inertial effects become more important. A possible interpretation 
could be that the decrease of T(X) for X < 1 is a signature of the adiabatic expansion of 
the ideal fluid when it climbs the potential barrier. Therefore it should be restricted to large 
C and low T, to avoid the thermalization with the bath. 

Finally, we present in Fig. fig:7 the results for the total force Ft, from eq. fl33|) . obtained 
both for the low and high shifting rates, as functions of the damping T. We present the 
results for three different choices of the parametrization, [i max from 2 to 4, so that they give 
also a picture of the convergence of our treatment in terms of the free modes of the system. 
Notice that from eqs. ( 1331) and ( 1351) . the same results may be scaled to get the total power 
pumped by the barrier, and are directly associated to the excess of mass, and of kinetic 
energy through the relationship (|36|) imposed by the steady state condition. The results 
shows a clear difference between the C = 0.2 case, with very little dependence of Ft on T, 
and the high rate shift, C = 2., with a very rapid decay of the force for increasing damping. 

V. CONCLUSIONS 

Colloidal particles when subjected to external driving forces exhibit may properties which 
are different from those of equilibrium systems. In the present paper we have described the 
effect of a barrier moving at constant velocity in a one dimensional colloidal fluid in the 
approximation that the solvent is unaffected by the barrier. In contrast with previous 
approaches which have considered only overdamped dynamics, we have studied the case 
where inertia plays a role. The two major effects of inertial terms are first to determine the 
appearance of a wake structure, completely absent in the DDF treatment and of an infinite 
set of characteristic lengths in the regions near the moving barrier; and second to produce not 
only a strongly structured density distribution near the barrier, and the associated current 
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density, but also the higher order moments of the velocity distribution distribution, which 
may be represented as a local temperature profile, very different from that of the thermalizing 
bath, and which shows interesting characteristics. It is also interesting that the method used 
here, based on the natural expansion of the distribution P(X, V, r) for free systems, gives 
an intuitive connection with the previous results based on the DDF treatment, i.e. using 
the density distribution p(X,r) = p o $ (X, r) as the only relevant field. That approach is 
recovered in as the limit of the simplest description of P(X, V, r) in terms of the first free 
mode (the only one with a purely diffusive dynamics, without an exponential decay time). 
The local balance of mass and force reproduce the DDF result of a Smoluchowski equation. 
To achieve the equivalent local balances for the energy, heat currents, etc., we have to 
enlarge the parametrization for P(X,V,t), to include exponential decaying modes, which 
represent the effects of the inertial dynamics of the particles. 

We have analyzed here only the simplest case, of one- dimensional spatial distributions 
in the dilute, ideal gas, limit. The equivalent results under other geometrical conditions, 
when the particles can bypass the moving barrier— ^ and including the effects of the particle 
interactions^ 1 ^, have been explored under the DDF assumptions, and it would be interesting 
to generalize them to the present approach. 

It is perhaps worth to comment that the wake region is not specific to the flow of particles 
with inertia. Very similar effects were also found in higher dimensions for overdamped, 
Brownian particles driven past colloids, which act in this case as the potential barriers^ 2 . As 
shown by Penna et al.— there exists a sum rule stating that the integral for the wake in any 
transversal plane to the direction of the drift vanishes, so that the depletion along the axis 
through the obstacle is exactly cancelled by the contribution from the lateral wings. Such a 
sum rule is valid in any dimension, but of course, in D = 1 implies that there is no wake at 
all. The presence of a wake structure in D = 1, in the inertial case, should correspond to a 
breaking of the sum rule for its transverse integral in D > 1. 

The present results perhaps are of relevance for microfluidic devices where colloidal par- 
ticles move along narrow channels in order to understand what external forces are needed 
to induce a drift in the presence of Brownian fluctuations. The hydrodynamics interaction, 
which have been neglected here, could not be important in one dimension due to the screen- 
ing effect, and we may expect that, at least at a qualitative level, the predictions made here 
could be accessible to experimental observation. 
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Caption List 
Caption Fig. 1 

Steady state scaled density profile, $oPO = PpO / Po, in the reference frame of the moving 
parabolic barrier. The barrier strength is k = 10, its width is 2 and moves at a relatively low 
velocity, C = 0.2, while the damping constant T takes on several values. Panel (b) shows 
the structure of the region within the barrier (—1 < X < 1), and the inertial wake left 
behind by the advancing barrier. The position X is relative to the barrier. Adimensional 
units (jlJJ-Q are used for all the quantities. 

Caption Fig. 2 

Scaled density profile, &o(X) = p(X)/p Q , at the front (a) and wake (b) regions in Fig.(l) 
is presented in reduced distances (X ± 1)CT, to take into account the natural decay length 
of the zeroth order mode. The front structure curves collapse for T > 1, while the wake 
region is reduced for increasing T. The position X is relative to the barrier. Adimensional 
units (H])-© are used for all the quantities. 

Caption Fig. 3 

Steady state density distribution, $o(A) = p(X)/p Q , induced by the same potential 
barrier as in Fig. |Vj but moving at a higher velocity, C = 2. The damping constants are 
T = 5 (full line), T = 4 (long dashed line), T = 3 (dot-dashed line), T = 2 (short dashed 
line). Panel (a) gives a general view of the high density structure at the advancing front. 
Panel (b) shows the structure of the depleted region within the barrier (—1 < X < 1), and 
the inertial wake leaved behind by the advancing barrier. The position X is relative to the 
barrier. Adimensional units (Jl])-© are used for all the quantities. 

Caption Fig. 4 

The structure of the relative density § (X) = p(X)/p , at the front (a) and wake (b) 
regions in Fig. (3) is presented in reduced distances (X ± 1)CT, to take into account the 
natural decay length of the zeroth order mode. The position X is relative to the barrier. 
Adimensional units (BJ-® are used for all the quantities. 

Caption Fig. 5 

Steady state temperature profile induced by a parabolic potential barrier shifted at rate 
C = 0.2. The position X is relative to the barrier and the vertical dotted lines are the 
barrier edges. The inset shows the structure of for T(X) « T Q . Adimensional units (SJ)- 
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are used for all the quantities. 
Caption Fig. 6 

Steady state temperature profile induced by a rapid drift, C = 2., of the parabolic 
potential barrier. The position X is relative to the barrier and the vertical dotted lines are 
the barrier edges. Adimensional units are used for all the quantities. 

Caption Fig. 7 

Steady state total force exerted by the moving barrier on the particles for the shifting 
rates C = 0.2 and C = 2. The full lines are obtained using the expansion (fl9|) in terms of 
the free modes up to fi ma x = 4, the dashed lines up to fi m ax = 2 and the dotted lines up to 
Umax = 2. Adimensional units eqs.dl])-© are used for all the quantities. 



21 



Fig. 1 
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Fig. 2 
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Fig. 4 
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Fig. 5 
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Fig. 7 
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